Effects of a nano-sized filler on the structure and dynamics of a 
simulated polymer melt and the relationship to ultra-thin films 



Francis W. Starr, Thomas B. Schr0der, and Sharon C. Glotzer 
Polymers Division and Center for Theoretical and Computational Materials Science, National 
Institute of Standards and Technology, Gaithersburg, Maryland 20899, USA 

(19 July 2000) 



Abstract 

We perform molecular dynamics simulations of an idealized polymer melt 
surrounding a nanoscopic filler particle to probe the effects of a filler on the 
local melt structure and dynamics. We show that the glass transition tem- 
perature T g of the melt can be shifted to either higher or lower temperatures 
by appropriately tuning the interactions between polymer and hller. A grad- 
ual change of the polymer dynamics approaching the filler surface causes the 
change in the glass transition. We also find that while the bulk structure of the 
polymers changes little, the polymers close to the surface tend to be elongated 
and flattened, independent of the type of interaction we study. Consequently, 
the dynamics appear strongly influenced by the interactions, while the melt 
structure is only altered by the geometric constraints imposed by the pres- 
ence of the filler. Our findings show a strong similarity to those obtained for 
ultra-thin polymer films (thickness ^ 100 nm) suggesting that both ultra-thin 
films and filled-polymer systems might be understood in the same context. 

PACS numbers: 61.25.Hq, 61.43.Fs, 61.82.Pv, 64.70.Pf, 66.10.Cb, 68.35.Ja 
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Significant enhancements in mechanical, rheological, dielectric, optical, and other prop- 
erties of polymer materials can be obtained by adding fillers such as carbon black, talc, silica, 
and other inexpensive materials [[I]]. Applications of filled polymers are diverse, ranging from 
automobile tires and bumpers to the rapidly expanding area of micro- and nano-electronic 
devices, which are continually challenged by the progressive miniaturization of devices [|T],3 . 
In particular, the growing ability to design customized nano-fillers of arbitrary shape and 
functionality provides an enormous variety of property modifications by introducing specific 
heterogeneity at the nanoscale However, detailed knowledge of the effects of fillers on 

a polymer melt at the molecular level is lacking due to the difficulty of directly probing the 
polymer structure and dynamics in the vicinity of the polymer-filler interface. In this regard, 
molecular simulations provide an ideal opportunity for direct insight into filled materials. 
Additionally, the understanding of ultra-thin polymer films, which also have many impor- 
tant technological applications (e.g. paints, lubricants, adhesives, and electronic packaging), 
is a topic of continuing discussion [p-TEfl; the present results provide a framework in which to 
interpret experiments on filled polymers, and also possibly polymer thin films, which report 
both increases and decreases of the glass transition temperature T g |T7[], depending on the 
details of the system studied [|l|18|- 



Our findings are based on extensive molecular dynamics simulations of a single nanoscopic 
filler particle surrounded by a dense polymer melt (Fig. [1]). We simulate 400 chains of 20 
monomers each (below the entanglement length), a total of 8000 monomers. The polymers 
are modeled as chains of monomers, which interact via a Lennard Jones (LJ) potential 
Vlj(?") = 4e((<j/r) 12 — (cr/r) 6 ). Additionally, bonded monomers are connected via a FENE 
anharmonic spring potential Vfene = -k(H%/2) m(l-(r/i? ) 2 ) 0i3- Several authors have 
studied the pure (unfilled) melt |23|,|24[] , and the system has been shown to be a good glass- 



former. This type of "coarse-grained" model is frequently used to study general trends of 
polymer systems, but does not provide information for a specific polymer. The filler particle 
is modeled as an icosahedron, not unlike primary particles of graphitized carbon black |2B| . 
We assign ideal force sites at the vertices, at 4 equidistant sites along each edge, and at 6 
symmetric sites on the interior of each face of the icosahedron, as shown in the right panel 
of Fig. [I| We tether a particle to each of these sites by a FENE spring, which maintains 



a relatively rigid structure but allows for thermalization of the filler. p5| . We consider a 
filler particle with an excluded volume interaction only, as well as one with excluded volume 
plus attractive interactions (expected for many fillers), to determine which properties are a 
result of the steric constraints imposed by the filler, and which properties are affected by 
polymer-filler attraction. We choose the same parameters for the interaction potential for 
all filler force sites. Periodic boundary conditions are used in all directions. 

Our model filler has several general features typical of a primary carbon black particle 
(a traditional filler) HI],|26||, as well as some newer nano-fillers (i) it has a size of order 

10 nm and (ii) it is highly faceted, but nearly spherical. The size of the facets is roughly 
equal to the end-to-end distance R e of the low molecular weight polymers comprising the 
surrounding melt (for Gaussian chains, R\ = 6-R 2 , the radius of gyration). Since we aim to 
study the changes induced by a single filler on the nanoscopic level, and the implications of 
those changes for bulk properties, we also consider a pure dense melt for comparison. 

We simulate the pure system at density p = 1.0 at temperatures ranging from T = 0.37 
- 1.0 [28]. We simulate the filled systems in the range T = 0.35 - 1.2. Equilibration times 
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range from 5 x 10 2 r* at the highest T to 2 x 10 4 r* (approximately 40 ns in Argon units) 
at the lowest T; we use the rRESPA multiple time step algorithm to improve simulation 
speed [P?,p5 . In order to compare the simulations of the filled system with the pure melt, 



we choose the box size so that the local density far from the filler deviates at most by 
0.2% from the density of the pure melt; such a density difference would cause a change 
in T g in this model less than that reported in Fig. |] ||23|| . For attractive monomer-filler 
interactions, a box size L = 20.4 (all lengths in units of <7 mm ) satisfies this constraint at all 
T. In the non-attractive case, the characteristic first neighbor distance between the filler 
sites and monomers is T dependent due to the lack of a unique minimum in the polymer- 
filler interactions. As a result, at each T a different L is required to achieve the correct p at 
large distance from the filler. The box sizes range from L = 20.49 at T = 1.0, to L = 20.6 
at T = 0.4. 

To quantify the effect of the filler on T g — one of the most important processing parame- 
ters — and on dynamic properties, we first calculate the relaxation time r of the intermediate 
scattering function 

F (q, t ) = ( J2 e-^V~ r ^) , (1) 

\j,k=i 



NS(q) 



where the sum is performed over all N monomers and we normalize by the structure factor 
S(q) such that F(q, 0) = 1. We define the value of r by F(q, r) = 0.2. Relative to the pure 
system, we find that r is larger at each T for the attractive system (Fig. |^). Conversely, r 
appears slightly smaller at low T in the non-attractive system, but is nearly indistinguishable 
from the results for the pure system. The lines are fits to the Vogel-Fulcher-Tammann (VFT) 



form 32 



r ~ e A ' {T - T ^ (2) 



where T is typically quite close to the experimental T„ value [32]; hence changes in T g are 



reflected in T . Consistent with the changes in r relative to the pure melt, we find that T 
increases in the system with attractive interactions, but clearly decreases in the system with 
only an excluded volume interaction. Thus, the effect of the steric hindrance introduced by 
the filler particle decreases r(T) and T g , in spite of the fact that monomers have a reduced 
number of directions in which to move, and hence degrees of freedom that aid in the loss 
of correlations. The fact that T g shifts in opposite directions for attractive versus purely 
excluded volume interactions demonstrates the importance of the surface interactions. 

To elucidate how the local dynamics of the monomers are influenced by the filler, we 
examine the relaxation of the self (incoherent) part F se \f(q, t) of F(q, t) as a function of the 
monomer distance from the filler. Monomers typically form layers near a surface |l lj]; we 
find well-defined monomer layers surrounding the filler, as seen in the density profile of 
Fig. ||. Hence, we split F se if(g, t) into the contribution from each separate layer. Specifi- 
cally, we calculate F^ r (q,t) using the monomers located in each layer at t — such that 
F se ii(q,t) = 1/N Xaayers Ni a y ev F^ T (q, t), where iVi ayer is the number of monomers in a given 
layer. We show F s ^ I (q ,t), as well as F se n(q ,t) for one temperature in Fig. ||. In the 
attractive system, the relaxation of the layers closest to the filler are slowest, consistent 
with the system dynamics being slowed by the attraction to the filler. Conversely, for the 
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non-attractive system, we find that the relaxation of inner layer monomers is significantly 
enhanced compared to the bulk, consistent with the observed enhancement of the system 
dynamics. Preliminary results support the possibility that faster dynamics may also occur 
with attractive interactions, provided that the polymer-filler attraction is weaker than that 
of polymer-polymer interactions. The altered dynamics persist for a distance slightly less 
than 2R g from the surface. Our results demonstrate that interactions play a key role in 
controlling T g and the local dynamics of filled polymers. We expect the role of interactions 
to be largely the same when many filler particles are present in the melt, but there will be 
additional effects on dynamic properties due to the more complex geometrical constraints. 

We next turn our attention to any structural effect the filler has on the melt. Quantities 
such as the pair distribution function, average R g , average R e , and distribution of bond 
lengths and angles show no significant deviations from the pure system. However, by focusing 
on the dependence of R g (or R e ) on the distance d from the filler surface, we find a change 
in the overall polymer structure near the surface. In Fig. [|, we show R 2 , as well as the 
radial component from the filler center R^ 2 (approximately the component perpendicular 
to the filler surface) for both attractive and non-attractive polymer-filler interactions at one 
temperature. R 2 increases by about 50% on approaching the filler surface; at the same time 
Rg 12 decreases by slightly more than a factor of 2 for both attractive and non-attractive 
systems. The combination of these results indicates that the polymers become slightly 
elongated near the surface, and flatten significantly. Note that not all monomers belonging 
to a given "surface polymer" are located in the first surface layer, as depicted in Fig. [L| We 
also point out that the chains retain a Gaussian conformation near the filler surface J34j . We 
find that the range of the flattening effect roughly spans a distance R g from the surface, and 
the results depend only weakly on T. We performed an additional simulation with double the 
attraction strength between the filler and polymers and did not find any significant further 
effect on the chain structure. The independence of the chain structure on the interaction 
suggests that the altered shape of the polymers is primarily due to geometric constraints 
of packing the chains close (d ^ R g ) to the surface. For significantly stronger interactions, 
alteration of the chain structure is expected on theoretical grounds |35 |. 

We next consider the implications that our results may have for studies of ultra-thin 
polymer films, where there is long standing debate on the role of interactions versus con- 
finement on T g shifts [|5T, |i~2"| , |i"3| , local melt dynamics |5j,|l4]-|T6| , and melt structure |5|-|TT 



Our simulations allow us to address the effects of interactions with a surface, without the 
additional complication of confinement effects present in thin films. It is largely agreed 
that ultra-thin films with strongly attractive substrates increase T g , while weak substrate 
interactions (or no substrate, as in freely standing films) lead to a downward shift of T g ; 
this is consistent with our results. This consistency is reasonable for fillers which have 
facets that are relatively smooth and large compared to R g ; for nanoscopic fillers, such as 
we study, it is surprising that a correspondence occurs even for R g close to the filler size. 
Not surprisingly, the magnitude of the shifts we observe depends on the relative quantities 
of polymer and filler; a greater filler concentration would have a more dramatic effect (as 



observed experimentally in refs. [pjH 20|l). Insofar as the magnitude of effects depends only 
on the ratio of the surface to bulk monomers, the thickness of the film is analogous to the 
inverse of the concentration of the filler. This is consistent with the experimental observation 
that T g shifts are more pronounced as film thickness decreases. Recently there have been 
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several experiments on segmental motion in both freely-standing and supported ultra-thin 
films I), 14, 15 . The observed segmental dynamics is consistent with a decreased T g find 



in calorimetric measurements [^,0,|T^]. At this time, it is not clear whether a model with 
layers of different mobility is applicable to understand T g shifts of thin films ||; however, 
the parallel behavior we observe between the thin films and our simulations of a filled melt 
support this viewpoint. Finally, the elongation and flattening of polymers we observe near 
the filler has been observed in thin-film simulations |5|-|TT] as well as recent experiments 
H|7|]; the range of the effect found in Ref. || is quantitatively consistent with our results, 
which show the effect only for a range of roughly R g , while the results of Ref. J7| observed 
flattening for film thicknesses < 6R g . We also found, as in Ref. |J, that the chains retain a 
Gaussian structure near the surface. Thus our findings demonstrate that confinement is not 
a necessary ingredient for the observed changes in the dynamics and structure of polymers 
near surfaces. While our results provide strong support for interpreting the results for filled 
melts and ultra-thin films in the same framework, it is obvious that much care must be used 
when analyzing specific systems. 
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FIGURES 



FIG. 1. [Figure excluded because of file size]. "Snapshot" of our simulation of the filled polymer 
melt. The bonds between nearest-neighbor monomers along a chain are drawn in various shades 
of gray for clarity. The left panel shows the filler with the surrounding melt; the right panel shows 
a few representative polymers that have monomers near the filler surface. 




FIG. 2. Temperature dependence of the relaxation time of the collective intermediate scattering 
function. For attractive polymer-filler interactions, r is increased at all T relative to the pure 
system, and correspondingly has a greater value of To (and thus T g ). Conversely, r is reduced and 
To is smaller for non-attractive (excluded volume only) interactions. The lines are a fit to the VFT 
form. The inset shows the same data plotted against reduced temperature Tq/(T — To) to show 
the quality of the VFT fit. For clarity of the inset, r of the pure system is multiplied by 2, and r 
of the filled non-attractive system is multiplied by 4. 
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FIG. 3. The self (incoherent) part of the intermediate scattering function for the average of all 
monomers (dotted line) and decomposed into layers (defined by the distance from the filler surface) 
for (a) attractive interactions and (b) non-attractive interactions at temperature T = 0.4. The 
inset of each figure shows the local density profile p(d/R g ) of monomers as a function of distance 
from the filler, normalized by R g of the polymers. We define the distance d from the filler surface as 
the difference between the radial position of a monomer r mon and the radius of the inscribed sphere 
of the icosahedral filler particle ri cos = ^(42 + 18\/5) 5 L, where L is the length of an edge of the 
icosahedron. The monomers order in well-defined layers surrounding the filler; we use the minima 
in p(r) to define the boundary between layers. The main part of each figure shows F se \{(qo,t) for 
each of the layers. We calculate -F se if(<70) t) — given by Eq. |l| with the sum restricted to only one 
index — for monomers located in each layer at t = 0. In (a), we see that the relaxation near the 
filler surface is slowed by roughly 2 orders of magnitude. This is consistent with the dynamics 
being slower in the filled attractive system (and hence an increased T g ). In contrast, in (b) the 
relaxation of i*seif(9o>*) is enhanced near the surface. This is consistent with the dynamics being 
faster in the filled non-attractive system (and hence a decreased Tg). The relaxation time of the 
outer most layer in both cases nearly coincides with the relaxation time of the pure system. 
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FIG. 4. Radius of gyration R g of the polymer chains as a function of distance d/R g of the center 
of mass of a chain from the filler surface for T = 0.4. We also resolve the component perpendicular 
to the surface, which we label by R^. We show results for (a) attractive and (b) non-attractive 
interactions. The increase of R g as the surface of the filler is approached indicates that the chains 
become increasingly elongated. The decrease of R^ shows that the chains also become "flattened" 
in the radial direction - roughly perpendicular to the filler surface. The effect appears largely 
independent of the temperature and numerical values of the potential parameters. 
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